Thermopower of atomic-size metallic contacts 
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The thermopower and conductance of atomic-size metallic contacts have been simultaneously 
measured using a mechanically controllable break junction. For contacts approaching atomic dimen- 
sions, abrupt steps in the thermopower are observed which coincide with jumps in the conductance. 
The measured thermopower for a large number of atomic size contacts is randomly distributed 
around the value for large contacts and can be either positive or negative in sign. However, it is 
suppressed at the quantum value of the conductance Go = 2e 2 /h. We derive an expression that 
describes these results in terms of quantum interference of electrons backscattered in the banks. 

PACS numbers: 72.15.Jf, 72.10.Fk, 73.23.Ad, 73.40. Jn 

In recent years stable metallic contacts consisting of a single atom have become experimentally accessible Q. The 
interesting interplay between quantization of the electron modes and the atomic structure of the contacts has resulted 
in intensive research in this field. Information obtained from experiments on atomic-size metallic contacts has mainly 
been limited to the measurement of the conductance. Two exceptions stand out: the simultaneous measurements of 
force and conductance by Rubio et al. Q , which prove that the conductance steps produced by contact elongation are 
due to atomic rearrangements, and the measurement of the subgap structure in atomic size superconducting aluminum 
contacts which characterize the conduction modes, by Scheer et al. fl. In this paper we present measurements of the 
thermopower in atomic-size metallic contacts. 

The thermopower S is the constant of proportionality between an applied temperature difference A9 and the 
induced voltage, Vt p — SA8. The relationship between the thermopower and the electrical conductance G is given in 
the linear-response approximation by 
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with /i the chemical potential. One can view the thermopower as a measure for the difference in conductance between 
electron and hole quasiparticle excitations, or as the energy dependence of the conductance. We will argue that 
the dominant contribution to the thermopower in atomic size contacts comes from quantum interference terms as a 
result of backscattering of electrons on defects near the contact. We have studied the thermopower of atomic-size gold 
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FIG. 1. Schematic diagram of the modified MCB configuration, used for the simultaneous measurement of conductance and 
thermopower. 
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FIG. 2. Typical conductance G and thermopower S versus piezovoltage Vp. The vertical gray lines indicate the corresponding 
steps in the conductance and thermopower. For two plateaus the conductance scale has been expanded 10 times and offset in 
order to show the small anomalies in G. 



contacts using a mechanically controllable break junction (MCB) g. A schematic diagram of the sample configuration 
is shown in Fig.[l]. By bending the phosphor bronze substrate, the 100 /itm gold wire breaks at the notch, allowing 
atomic-size contacts to be adjusted. This gold wire is attached to long thin (25 /im) gold wires at both ends. They 
connect the notched wire to the current and voltage leads, anchored at the bath temperature, hence forming an open 
gold loop. The central gold wire is tightly wound and varnished on each side of the contact around a calibrated 5 kil 
R11O2 thermometer and a 500 Q RuC>2 heater. Using one heater, a temperature gradient over the contact can be 
applied. The glass plates and thin gold wires serve as thermal resistances to the substrate and bath temperatures, 
respectively. The sample is placed in a regular MCB setup in an evacuated can immersed in liquid helium. The 
conductance and thermopower were measured in three steps: The voltage over the contact was measured with a 
nanovoltmeter at — 100 nA, +100 nA and at zero dc current bias, while maintaining a constant temperature gradient 
over the contact by applying about 2 mW heating power to a heater on one side of the constriction. The conductance 
is then obtained from the voltage difference for the two current polarities, and S is obtained from the voltage at 
zero bias current. Each cycle takes about 4 s and is continuously repeated as we slowly sweep the piezovoltage up 
in order to decrease the contact size. A curve (Fig. ||) was generally taken from 10 Go to tunneling in 30 min (Go is 
the quantum conductance unit, 2e 2 /h). Every few traces, it was necessary to readjust the contact manually, which 
inevitably leads to the contact being pushed completely together. Low pass RC filters (10 Hz) were mounted in the 
circuit near the sample to prevent rectification of ac disturbances by the asymmetry of the voltage dependence of the 
conductance. Extensive measurements have been performed on two samples (referred to as samples 1 and 2). 

The primary limitation of the sample design is the thickness of the glass plates and thus the thermal insulation of 
the gold sample wire from the phosphor bronze substrate. The thickness is a trade-off between stability and thermal 
insulation. As a result of thermal currents flowing to the substrate, a thermal gradient is established in the sample 
wire between the thermometer and the contact. The measured temperature is hence not the actual temperature of the 
"hot" side of the constriction. We calibrate this thermal gradient by measuring the thermopower for large contacts, 
with resistances in the range 1 — 10 fi, as a function of heating power. First, we stress that the thermal resistance 
of the contact is orders of magnitude larger than that of the wire on either side, therefore the temperature difference 
over the contact can be taken to be independent of the contact size. This is corroborated by experiment, which shows 
that the mean value of the thermopower over the contact as a function of contact diameter in the range 0.1 - 100 SI 
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FIG. 3. Density plot of thermopower against conductance: The thermopower axis is divided into 0.125 fiV /K, and the 
conductance into O.lGo sections, creating 0.125 /iV/K by O.lGo bins. The occurrence of a particular combination of conductance 
and S is added to its corresponding bin, and the sum of 220 curves for the two samples is plotted in gray scale above. Black 
represents no data points and white more than 100. We note that above 10 Go fewer data have been taken. 



remains constant within an accuracy of 1%. We take advantage of the fact that for conventional point contacts the 
phonon drag contribution to the thermopower becomes negligible j^] . Since the contact is part of a uniform gold loop, 
the measured themopower corresponds to the phonon drag contribution of the leads only. The side that is not heated 
remains equal to the substrate temperature and we assume that the actual temperature difference over the contact 
is a fixed fraction, a, of the measured temperature difference. We then determine this fraction by comparing the 
measured large-contact thermopower as a function of temperature with literature values for the bulk thermopower of 
pure gold pJj}, which is nearly linear between 10 and 25 K, with a slope of ^0.05 /zV/K 2 . For the two gold samples 
discussed below, the model provides a good description when the fraction a is taken as 0.4 and 0.5, respectively. We 
estimate an error of about 20% for the temperature difference obtained. 

To have a reasonable signal level we need to apply a temperature difference of several kelvins. In the case of 
sample 1 the temperature difference A8 = 4K, and the average of the temperatures on both sides of the contact is 
6> au = 11.5K. For sample 2, A6 = 6K and av = l2K. As it is not obvious that this AO can be regarded as a small 
(linear) perturbation, we will take the full dependence of S on Ad into account in the analysis below. For all the 
results presented below, the bulk thermopower of the leads has been subtracted. 

While breaking the contact by increasing the piezovoltage, the usual plateaus in the conductance are observed 
,|JUJ^]. When heating one side of the contact we observe steps in the thermopower which occur simultaneously 
with conductance jumps from one plateau to the next. Each measured curve produces a different conductance and 
thermopower trace and a typical example is shown in Fig.||. Even tiny jumps or changes of slope of the conductance 
can be accompanied by large steps in the thermopower. On a conductance plateau, even though the conductance 
hardly changes, smooth variations in the thermopower are usually observed. Note that the thermopower of the 
contact can have both a positive or negative sign. When we do not heat, or heat both sides of the contact to the same 
temperature, we observe no thermopower voltage within the noise level of 300 nV peak to peak. 

In order to obtain statistical information about a possible correlation between the measured thermopower and 
conductance values, a density plot was constructed from the combined data of the 72 and 148 individual curves 
from sample 1 and 2, respectively (Fig.|^). The conductance axis was divided into 10 partitions per Go and the 
thermopower axis in partitions of 0.125 /iV/K. Then, the number of data points falling in each range of conductance 
and thermopower was counted and the results are represented in gray scale in Fig.|3j. In this figure we observe an 
increase in the spread of the thermopower values with decreasing contact size, with both positive and negative sign. 

A conductance histogram of the 220 curves is, within the statistical accuracy, in agreement with other gold his- 
tograms taken at low temperatures (e.g. Ref. 0|). Although the data are not presented here, similar results for the 
thermopower have also been observed in silver and copper samples, albeit for a more limited number of curves. 
The thermopower has a random value and sign, and seems to be much more sensitive to small changes in the 
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atomic geometry of the contact than the conductance. This is not expected from the simple adiabatic models for 
point contacts, which only predict a positive sign Jl^Jll[] . We propose an interpretation of this behavior in terms 
of coherent backscattering of the electrons near the contact: As a result of the interference of waves with different 
path length, the transmission of the contact will show fluctuations as a function of energy [|l2[|l3f |. Each atomic 
rearrangement at a conductance step will alter the interference paths of the backscattered electrons by a significant 
fraction of , and hence change the energy dependence of the transmission in an unpredictable way causing each step 
in the conductance to result in an unpredictable jump in the thermopower. Along a plateau, the contact gradually 
changes position with respect to the scattering centers nearby and a gradual change in the interference pattern occurs. 

We now derive an expression for the ther mopow er based on these concepts. The thermopower in quantum point 
contacts can be written as S = —L/G, with 
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where (Tr tt^) is the sum of the transmission probabilities, f(9, E) is the temperature- and energy-dependent Fermi 
function, h is Plank's constant, and e = |e| is the electron charge. The thermopower is characterized by the energy 
dependence of the transmission probabilities. We have approached the problem using the same method as presented 
in Ref. |14| for the derivation of the conductance fluctuations in atomic size contacts. The point contact is taken as a 
ballistic central constriction with diffusive regions on both sides. The ballistic section (using the Landauer-Biittiker 
formalism) is characterized by N conductance modes, each with a transmission probability T n and a contribution 
T n Go to the conductance. After transmission through the contact, within the dephasing time r^, the electron scatters 
elastically in the diffusive region and has a finite probability amplitude, a, to return to the contact. When the diameter 
of the contact is small compared to the elastic scattering length l e , the return probability, |a| 2 , is small and we need 
only consider lowest-order processes. To lowest order in a we can write the transmission of the three sections combined 
as 

Tr tt* = ^T n [l + 2Re {r n a hn + r' n a r<n )}. (1) 

Here, r n and r' n arc the reflection coefficients in the transfer matrix describing the central ballistic section of the 
contact when coming from the left and right, respectively, with \r n \ 2 = \r' n \ 2 = 1 — T n . a^ n and a r ^ n are the return 
amplitudes for mode n from the left and right diffusive regions, respectively. The latter are sums over all possible paths 
of length /, containing phase factors e l ( E - E F) l / v F h _ The second term in Eq. (|l|) describes the interference of the directly 
transmitted wave with the fraction that, after transmission, is first backscattered to the contact and subsequently 
reflected at the contact. Assuming the dominant energy dependence is in the phase factors, the integration over E 
in the expression for L can be performed. We consider the square of L, averaged over an ensemble of scattering 
configurations, 

with z = irkBOr/h and r = l/vp. Here we have assumed that a\ n and a r n are uncorrelated and have the same 
average return probability (|a| 2 ), independent of the mode number n. 

For (|a(r)| 2 ) we substitute the semiclassical probability to return to the contact into any of the TV modes after a 
time r, (|a(r)| 2 ) = vp /{\/12Trkp(DT) 3 / 2 ], with D = vfIc/3 the diffusion constant. The integral in Eq. (0) can be 
performed numerically. It only weakly depends on the ratio A9/9. For the standard deviation of the thermopower 
a s = y/(S 2 } - (S) 2 = ^(L 2 )/(G 2 ) this finally results in 
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a s = - , — [ ) ^ — "1 (3) 



Here, c is a numerical constant which equals 5.658 in the limit A9/9 — * 0, and increases by about 5% for A9/9 = 0.5. 
We have also introduced a factor (1 — cos 7) to account for the finite geometrical opening angle of the contact, where 
the limit 7 = 90° corresponds to an opening in an infinitely thin insulating layer between two metallic half spaces. 
Note that as in Eq. (||) is equal to zero when all T n are equal to either or 1 . 

In Fig.|] we plot the standard deviation of the thermopower, determined from the experimental data by sorting 
all data points as a function of G from the combined set of the 220 curves and averaging over 1000 consecutive data 
points. We compare these data to the theoretical curve calculated from the above expression for the case where the 
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FIG. 4. Standard deviation of the thermopower against conductance. Solid squares represent the measured data and the 
solid line is the theoretical curve assuming the conductance modes open one by one. 



modes contributing to the conductance open one by one. That is, for all conductances with TV modes contributing 
to the conductance, only one (i.e., T/v) differs from 1 and as oc T/vVl — Tn/(N — 1 + T/v). The observed deep 
minimum at G = Go suggests that this conductance is dominantly carried by a single mode. This is in agreement 
with measurements of the shot noise on atomic size gold contacts [^5j, and the observed suppression of conductance 
fluctuations in Ref. fji|j . For G > Gq, the limited statistics in combination with the property that the effect in the 
thermopower scales inversely with conductance prevent the definite identification of minima near quantized values. 

From the amplitude of the curve we obtain an estimate for the elastic mean free path of l e = 5 ± 1 nm, using 
reasonable values for the opening angle of the contact of 35°-50° jl6| . All data points should be on or above the full 
curve in Fig. ^, since contributions by more conductance modes can only increase the variance in the thermopower. 
Therefore, l e cannot be much smaller than 4 nm. For a much larger value of l e many modes would have to contribute 
to the conductance, in which case we would not expect to find any minima at quantized values. 

Apart from the thermopower effects described here, for a quantum point contact positive peaks in the thermopower, 
centered at conductance values (n + ^)Go, n = 0,1,2,..., were predicted due to the structure of the electron density of 
states • This effect has indeed been observed in two-dimensional electron gas devices 1 17 , but is much smaller 



than the fluctuations we observe here, and therefore cannot be resolved in the mean value (S) for our metallic point 
contacts. 

The mechanism we present to explain the thermopower is the same as the one proposed for the voltage dependence 
of the conductance Indeed, when we plot as and y/ {(dG / dV) 2 ) / G for gold, the data show very similar behavior. 
The energy scales with which both measurements have been performed are so different (6 K temperature difference 
versus 20 mV amplitude) that comparison between the parameters obtained by both methods is a test for the validity 
of the theoretical derivation. The mean free path value obtained from the conductance fluctuations is 5 nm, in good 
agreement with the estimate obtained here. Hence, the fact that both works are not only in good qualitative but also 
in good quantitative agreement is strong support for the model used. 

This work was part of the research program of the "Stichting FOM" , which was financially supported by NWO. We 
have greatly profited from many discussions with C. Urbina, D. Esteve, M.H. Devoret, and we thank L.J. de Jongh 
for his stimulating support. 
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